@hfhsieh

Since it's only during initialization that we do not need to update the magnetic energy in MHD_SameInterfaceB(), I think we can always perform the update, as it should not impact performance.

Agreed. I’ll let you know once I’ve updated the code.

Should we update the wiki documentation to note this behavior under MINMOD_MAX_ITER, or consider automatically enabling OPT__SAME_INTERFACE_B when MINMOD_MAX_ITER > 0 in MHD simulations?

I will update the wiki document.
@hyschive Would you like to automatically enabling OPT__SAME_INTERFACE_B when MINMOD_MAX_ITER > 0 in MHD simulations?

I will update the wiki document. @hyschive Would you like to automatically enabling OPT__SAME_INTERFACE_B when MINMOD_MAX_ITER > 0 in MHD simulations?

Agree. How about declaring OPT__SAME_INTERFACE_B as int instead of bool, setting its value to -1 in all Input__Parameter, and resetting it to either 0 or 1 (if its value is still -1) in Init_ResetParameter.cpp depending on the value of MINMOD_MAX_ITER?

You can also define a corresponding typedef in Typedef.h for OPT__SAME_INTERFACE_B.

Also, remember to raise a warning message in Aux_Check_Parameter.cpp when enabling MINMOD_MAX_ITER for MHD.

@hfhsieh I’ve updated the PR. Please take a look when you have time—thank you in advance!

@hfhsieh Thanks for the review. I have addressed the comments.

@ChunYen-Chen Could you resolve the conflicts caused by the recent wiki updates? Thanks.

@ChunYen-Chen Is it ready for another review?

@hyschive There is one last comment that needs to be resolved. I will let you know once I have finished.

@hyschive Thanks for the review. The PR has been updated and is ready for the next round of review.